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O ABSTRACT 
CN 

O Aims. Studying the effects of an inhomogeneous interstellar medium (ISM) on the strength and the shape of the Lyman 

alpha (Lya) line in starburst galaxies. 

Methods. Using our 3D Monte Carlo Lya radiation transfer code, we study the radiative transfer of Lya, UV and optical 
continuum photons in homogeneous and clumpy shells of neutral hydrogen and dust surrounding a central source. Our 
simulations predict the Lya and continuum escape fraction, the Lya equivalent width EW(Lya), and the Lya line 
profile, and their dependence on the geometry of the gas distribution and the main input physical parameters. 
Results. The ISM dumpiness is found to have a strong impact on the Lya line radiative transfer, entailing a strong 
dependence of the emergent features of the Lya line (escape fraction, EW(Lya)) on the ISM morphology. Although a 
clumpy and dusty ISM appears more transparent to radiation (both line and continuum) compared to an equivalent 
homogeneous ISM of equal dust optical depth, we find that the Lya photons are, in general, still more attenuated 
• than UV continuum radiation. As a consequence, the observed equivalent width of the Lya line (i5IVobs(Lya)) is lower 

i^H than the intrinsic one (i5IVint(Lya)) for nearly all clumpy ISM configurations considered. There are, however, special 

conditions under which Lya photons escape more easily than the continuum, resulting in an enhanced iJWobs(Lya). 
JL The requirement for this to happen is that the ISM is almost static (galactic outflows < 200 km s"'^), extremely clumpy 

^ (with density contrasts > lO'^ in HI between clumps and the interclump medium), and very dusty (E(B-V) > 0.30). 

When these conditions are fulfilled the emergent Lya line profile generally shows no velocity shift and little asymmetry. 
^ Otherwise the Lya line profile is very similar to that expected for homogeneous media. 

Conclusions. Given the asymmetry and velocity shifts generally observed in star-forming galaxies with Lya emission, 
we therefore conclude that clumping is unlikely to significantly enhance their relative Lya/UV transmission. 
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1. Introduction strong Lya emission line. Unfortunately, the first attempts 

to detect high-redshift galaxies in Lya gave quite meager re- 

^ Being the intrinsically brightest spectral signature of re- suits. The observed Lya fluxes appeared fainter than those 

^ mote young galaxies (Partridge & Peebles 1967; Schaerer predicted and only few Lya emitters (LAEs) had been de- 

^ 2003), and possessing a rest wavelength of 1216 A (making tected until the late 1990s (cf. Djorgovski & Thompson 

IJ it accessible for optical/near-IR ground-based telescopes 1992, Pritchet 1994). This lack of Lya emission has nev- 

.5^ for redshifts z > 2), the Lyman alpha (Lya) line has be- ertheless triggered several studies which have enabled to 

^ come the most powerful emission-line probe of the distant highlight the high complexity of the resonant Lya line ra- 

5— ( young universe. The potential of the Lya emission line diative transfer in starburst galaxies (Meier & Terlevich 

for detection and redshift confirmation of distant galax- 1981; Neufeld 1990; Chariot & Fall 1993; Kunth et al. 1998; 

ies, derivation of star formation rate (SFRs), as well as Tenorio-Tagle et al. 1999; Mas-Hesse et al. 2003; Ostlin 

probe of the ionization state of the intergalactic medium et al. 2009). While the faint measured Lya fluxes were orig- 

(IGM) (Malhotra & Rhoads 2004; Kashikawa et al. 2006) inally attributed to the dust attenuation (Pritchet 1994), 

and the reionisation epoch (Fan et al. 2002; Santos 2004) it has turned out that many physical effects could strongly 

is enormous, but necessarily reUes on a good astrophysical modify or supress the Lya line within galaxies (metallicity, 

understanding of the processes that regulates the emergent neutral hydrogen kinematics, geometry of the interstellar 

Lya emission from a galaxy. medium (ISM)). It is only during this last decade, and the 

The importance of the Lya line in the cosmological con- development of deep and wide surveys, that many Lya- 

text was first proposed by Partridge & Peebles (1967) who emitting galaxies have been detected (Hu et al. 1998, 2004; 

suggested that young high-z galaxies, undergoing their first Cowie & Hu 1998; Kudritzki ct al. 2000; Rhoads et al. 2000; 

star-forming event, should be detectable thanks to their Taniguchi et al. 2003, 2005; Shimasaku et al. 2006; Gronwall 
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ct al. 2007; Nilsson et al. 2007; Guaita et al. 2010; Ouchi 
et al. 2003, 2008, 2010). 

Because of the factors which contribute to the Lya ra- 
diative transfer, the Lya line features (line profile, equiva- 
lent width EW(Lya), offset from other emission/absorption 
lines) encode much information on the properties of individ- 
ual galaxies: gas kinematics, gas geometry as well as stellar 
population. For instance, the detection of unusually strong 
Lya line in the spectra of high-z galaxies could indicate 
the presence of population III stars within them (Schaerer 
2003), whereas the asymmetry of the line profiles would 
suggest the presence of strong galactic outflows (Kunth 
et al. 1998). The derivation of this precious information 
requires however an accurate interpretation of the Lya line 
features, which implies as well a complete understanding of 
the Lya radiative transfer in the ISM of galaxies. This is 
one of the aims of this paper. 

Due to the importance of the Lya line for cosmology, 
several studies have attempted to understand the physical 
process governing the escape of Lya photons from galaxies. 
Among the parameters which influence the visibility of the 
Lya line, dust content, neutral gas kinematics and the ge- 
ometry of the neutral gas seem to play the most important 
roles. Dust was originally invoked to explain the absence 
or the faint Lya emission from galaxies at large redshift 
(e.g. Meier & Terlevich 1981). However, Giavalisco et al. 
(1996) studied a local sample of star-forming galaxies ob- 
served with the lUE space telescope and found no clear cor- 
relation between Lya/H^ or EW(Lya) and the reddening 
E{B—V). Other studies have also led to a lack of correlation 
between the dust attenuation and the strength of the Lya 
line, suggesting that other parameters govern the escape of 
Lya photons (Kunth et al. 1994; Thuan et al. 1997; Atek 
et al. 2008). Among them, the role of the neutral gas kine- 
matics was revealed in the 1990s by Kunth et al. (1994) 
and Lequeux et al. (1995). For 8 local galaxies observed 
with the Goddard High Resolution Spectrograph (GHRS), 
Kunth et al. (1998) found that when Lya line appeared in 
emission, there was a systematic blueshift of low ionisation 
states (LIS) metal absorption lines with respect to Lya, in- 
dicative of outflows in the neutral medium. Furthermore, 
the shape of the Lya line profiles proved to be asymmetric. 
Galaxies showing Lya in absorption showed significantly 
smaller relative shifts of LIS lines and Lya. This result 
clearly shows that the Lya escape fraction and line shape 
are strongly affected by the kinematical configuration in 
the ISM. Phenomenologically it is easy to understand that 
an outfiow in the neutral ISM would promote the escape of 
Lya photons and create asymmetric line profiles, since the 
motion Doppler shifts the line out of resonance and more 
so for the red side of the line. Finally, several studies of the 
resonant Lya transfer have emphasized the importance of 
the ISM dumpiness on the escape of the Lya (Neufeld 1991; 
Giavalisco et al. 1996). In particular, Neufeld (1991) showed 
that it could be possible to observe an emergent EW(Lya) 
higher than the intrinsic one in a dusty and clumpy ISM. As 
the dumpiness of the ISM is well established in our galaxy 
(Stutzki & Guesten 1990; Marschcr et al. 1993), this pa- 
rameter must therefore be taken into account in the study 
of the Lya radiative transfer. 

With the increased number of Lya radiative trans- 
fer codes developped recently (Ahn et al. 2001, 2002; 
Cantalupo et al. 2005; Verhamme et al. 2006; Pierleoni et al. 
2007; Laursen et al. 2009a; Forero-Romero et al. 2011), the 



transfer of Lya photons has intensively been investigated 
in the framework of galaxy simulations. In particular, such 
simulations allow us to compare the observed Lya line prop- 
erties of individual galaxies, both nearby and distant ones 
(Ahn et al. 2003; Verhamme ct al. 2008; Atek et al. 2009). 
However, although most studies have treated the Lya radia- 
tive transfer in either static or expanding media, the main 
eff'ects of a multiphase ISM on the Lya radiative transfer 
has been the object of few numerical studies (Haiman & 
Spaans 1999; Richhng 2003; Hansen & Oh 2006; Laursen 
et al. 2012). The aim of this present paper is to carry out 
a detailed study of both the Lya and the UV continuum 
radiative transfer in a large range of dusty, moving, homo- 
geneous and clumpy ISMs. This will allow us to examine 
in detail the effects of the ISM dumpiness on the features 
of the Lya line (Lya escape fraction, EW(Lya), Lya line 
profiles) . 

One of the main motivations of our study is also to un- 
derstand the anomalous strong EW(Lya) revealed by sev- 
eral observations of LAEs at high-z (Kudritzki et al. 2000; 
Malhotra & Rhoads 2002; Rhoads et al. 2003; Shimasaku 
et al. 2006; Kashikawa et al. 2012). While normal stellar 
population models predict a maximum value of ~ 240 A 
for the intrinsic EW(Lya) within starburst galaxies (i.e. as- 
suming population I/II stars. Chariot & Fall 1993; Schaerer 
2003), it is not rare to observe higher EW(Lya) from high- 
redshift sources. Several physical possibilities have already 
been investigated to explain these high EW(Lya), such 
as the presence of either population III stars or Active 
Galactic Nuclei (AGNs) in the host galaxies. But none of 
them prove to be consistent with the observations (Dawson 
et al. 2004; Wang et al. 2004; Gawiser et al. 2006). Another 
possibility is that the high EW(Lya) values found are due 
to the combined effect of IGM absorption (lowering the 
continuum on the blue side of Lya at high z) and obser- 
vational errors biasing the average EW(Lya) to higher val- 
ues (Hayes & Ostlin 2006). The most popular explanation 
seems, however, to be the relative boost of Lya photons 
result in a clumpy ISM as originally suggested by Neufeld 
(1991). In this scenario, Lya and UV continuum photons 
propagate in a clumpy ISM, where all neutral hydrogen 
and dust are mixed together in clumps. While Lya pho- 
tons would scatter off of the surface of clumps, having their 
journey confined to the dustless interclump medium, the 
UV continuum photons would penetrate into the clumps 
and would suffer greater extinction. Such a scenario would 
thus produce larger EW(Lya) than the intrinsic ones, al- 
lowing to explain the anomalously high EW(Lya) observed 
in some high-z galaxies. Other studies have also invoked a 
higher transmission of Lya photons than for the UV con- 
tinuum, to explain observations of some low redshift Lya 
emitters (Scarlata et al. 2009), to understand the overall 
SED of LAEs at z ~ 4 (Finkelstdn et al. 2008, 2009), and 
to reproduce the Lya and UV luminosity function of distant 
galaxies (Dayal et al. 2009; Forero-Romero et al. 2011). 

In this paper, we investigate the Neufeld scenario fur- 
ther and examine the physical conditions under which a 
clumpy ISM could produce a boost of the Lya line relative 
to the continuum. 

The remainder of this paper is structured as follows. In 
Section 2 we outline a description of our numerical model, 

presenting the features of the clumpy media and our as- 
sumptions. In Sections 3 and 4 are presented our results. 
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The Lya radiative transfer in homogeneoiis and clumpy me- 
dia is presented in Section 3, whereas the formation and the 
features of the emergent Lya hne profiles are described in 
Section 4. Section 5 is dedicated to the discussion of these 
results with, in particular, an application of our study to 
the Neufeld scenario. Finally, our main conclusions are sum- 
marised in Section 6. 



2. Method 

2.1. 3D radiation transfer code 

To study the Lya line and UV-optical continuum radiation 
transfer in clumpy geometries, we have used the latest ver- 
sion of the 3D Monte Carlo radiative transfer code MCLya 
of Verhamme et al. (2006) and Schaerer et al. (2011). To 
treat the radiation transfer at wavelengths other than Lya, 
we here also compute the continuum transfer at other wave- 
lengths assuming scattering and absorption by dust. For the 
present paper we are interested in three wavelengths, listed 
in Table 2: the Lya line (A = 1215.67 A) and its neighboring 
UV continuum, and the optical B and V bands. 

2.2. 3D geometries, model parameters, and model output 

Both for simplicity, and since spherically symmetric out- 
flows with a homogeneous Hi shell are able to reproduce a 
large variety of observed Lya line profiles in Lyman break 
galaxies and Lya emitters (Verhamme et al. 2008; Schaerer 
& Verhamme 2008; Dessauges-Zavadsky et al. 2010), the 
same geometry is used to study how a clumpy ISM struc- 
ture alters the Lya line and UV continuum. This clumpy 
geometry is also chosen since it has been shown to repro- 
duce observable continuum properties of starburst galaxies 
and the Calzetti attenuation law (Gordon et al. 1997; Witt 
& Gordon 2000; Vijh et al. 2003). Finally, this also allows 
us to make a detailed study into the continuity of the ex- 
tensive grid of radiation transfer models by Schaerer et al. 
(2011). 

In practice we adopt the following, simple shell geome- 
tries (see figure 1): a static or radially expanding, homoge- 
neous or clumpy shell of Hi and dust surrounding the source 
emitting both Lya line and continuum photons. Dust and 
gas (Hi) are assumed to be co-spatial in the shell. We as- 
sume a point-like central source. 



2.2.1. Input parameters 

The four physical and the two geometrical input parame- 
ters of our models, listed in Table 1, are the following. The 
radial expansion velocity Vexp, the Doppler parameter h of 
the Hi, the mean Hi column density A^hi; the mean dust 
absorption optical depth Tq". the clump volume filling fac- 
tor FF, and the density contrast nic/«c between the inter- 
clump and clumpy medium. Each parcel of the shell (clump 
or interclump) exhibits the same radial velocity Vcxp- The 
Doppler parameter h — w^j^ -I- v^^^^ reflects the random 
(thermal -f turbulent) motions of the Hi. The clumpy (in- 
homogeneous) medium is defined by the volume filling fac- 
tor FF of clumps, by their density nc, and by the density 
contrast uiq / nc between clumps and interclumps of lower 
density nic- The mean Hi column density A^hi is thus re- 



lated to the (inter)clump density, FF, and the thickness of 

the shell L by: 



iVm = (FFnc + (1 - FF)nic)i. 
Similarly one has 

Ta = (1 - a)ad{ ){ — )Nm, 

nid Mh 



(1) 



(2) 



where a is the dust albedo, Od the total dust cross section 
(scattering -|- absorption), mjy the proton mass, the 
dust grain mass, and {Md/Mn) is the dust-to-gas ratio. In 
the present paper, the dust optical depth fa — the single 
parameter used to vary the dust content — is derived as- 
suming a dust grain size of 2xl0~^ cm and a mass rud = 
3x10"^'' g. The total dust optical depth is defined as 



Td 



(1 



(3) 



(4) 



and the dust particle density is 

m// Md , 
rid = nH — hr^ , 
rud Mh 

where nu stands for the clump or interclump density. We 
adopt the SMC dust properties (albedo a and phase func- 
tion g) listed in Table 2. These properties, together with the 
clumpy shell geometry also adopted here, have been shown 
to reproduce observable continuum properties of starburst 
galaxies (Gordon et al. 1997; Witt & Gordon 2000; Vijh 
et al. 2003). Although detailed model predictions depend to 
some extent on the dust properties, the main quantities of 
interest in this paper - the Lya and UV continuum escape 
fractions, and especially their relative values - should not 
strongly depend on the exact dust properties. We expect 
that other poorly known properties such as the geometry 
and velocity field, known to affect sensitively the transfer of 
Lya radiation, are more important than the detailed dust 
properties (Laursen et al. 2009b). For these reasons we have 
not considered changes of the dust properties, but focus on 
the effect of geometry and dumpiness in this paper. 

To construct clumpy structures with the desired input 
parameters in practice, we follow a similar approach as Witt 
& Gordon (2000). We construct a Cartesian grid of A^ = 
128'^ cells, within which the shell of thickness L is defined 
by an inner and outer radius, Rmin and i?max- Assuming a 
density nc for the high density regions (clumps) we then 
randomly choose a fraction FF of the cells localised in the 
shell (i.e. cells localised at a radius R such as i?min ^ R 
< ^max), which receive a high density nc- The remaining 
cells in the shell are set to low density njc- The physical 
cell size (or equivalently L) is then adjusted to reproduce 
the desired mean radial Hi column density A^hi, which is 
computed by drawing random lines of sight through the 
shell. Finally the dust content is varied by changing the 
dust-to-gas ratio [Md/Mn), yielding different values of the 
mean dust absorption optical depth Ta- 
in Table 3 we summarise the different values that we 
have explored for the six input parameters describing our 
models. For the present study we have adopted a filling 
factor FF = 0.23, as explained below (Sect. 2.2.2). The 
density contrast nic/nc has been varied from 1 (homoge- 
neous medium) to 0, reflecting the extreme case of an empty 
interclump medium. Models have been computed for static 
shells (vexp = 0) and expansion velocities up to Vexp = 600 
Then, a wide range of parameter space has been 



km 

considered, as listed in Table 3. 
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Homogeneous Clumpy Extremely clumpy 




Fig. 1. Representation of some 3D homogeneous and clumpy geometries studied in this paper. The star distribution is always 
localised in the center of the shell, whereas the dust and the Hi content are distributed around. The dust and Hi distribution 
can be "homogeneous" (left), "clumpy" (middle) or "extremely clumpy" (right). In a "clumpy" distribution, the clumps and 
the interclump medium receive, respectively, a high and a low densities of dust and HI. In the case of an "extremely clumpy" 
distribution, all the dust and the Hi content are distributed in clumps. 



Table 3. Range of values of the six input parameters (column 1-6) describing the homogeneous and clumpy shell models, and 
derived properties (cols. 7, 8). 



FF nic/nc 


Vcxp[km s ■^J b [km s 


iVm [cm-2] 


Ta 


CF mass spectrum 


0.23 [0, 1] 


0, 50, 100, 200, 250, 300, 400, 600 12.8, 20, 40 


[10^', 10^"] 


[0, 100] 


[0.7, 1] p(m) oc m-''-"" 



Table 1. Six input parameters (top) for and derived parameter 
(bottom) of the homogeneous and clumpy shell models 



Parameter 


Symbol 


Radial velocity 


Vexp 


Hi velocity dispersion 


6 


Mean Hi column density 




Mean dust absorption optical depth 




Clump volume filling factor 


FF 


Density contrast 


mc/nc 


Covering factor 


CF 



Table 2. Dust parameters (a and g) taken from Witt & Gordon 
(2000) and adopted for Lya line photons, and continuum pho- 
tons at UV and optical wavelengths (close to the B and V-band). 



Photons 


A (A) 


Td/rv 


a 


g 


Lya 
UV 


1215.67 
1235.0 


6.74 


0.460 


0.770 


B-band 


4350.0 


1.38 


0.495 


0.633 


V-band 


5550.0 


1.00 


0.490 


0.607 



2.2.2. Characterisation of clumpy structures 

Given a choice of the clump volume filling factor FF and the 
thickness of the shell (i.e. i?max — ^min), two other inter- 
esting quantities describing the inhomogeneous structure 
can be derived. First, the covering factor CF of the shell 
corresponding to the fraction of solid angle covered by the 
clumps as seen from the central source. Models with differ- 
ent covering factors are constructed by varying i?min- 



Another interesting quantity is the mass spectrum of the 
clumps. As clumps we consider, as Witt & Gordon (1996, 
2000), all cells directly connected with each other by at least 
one face. We then determine their mass spectrum, which 
approximately follows a power law p(m) cx m~", where m 
is the clump mass. Adopting FF = 0.23, we obtain a power 
law p{m) oc m~^ °^ as illustrated in Fig. 2. This mass spec- 
trum is consistent with observations of diffuse interstellar 
clouds showing a power law with a = 2 (Dickey & Garwood 
1989). The value FF=0.23 in our model is then the most 
appropriate value if we aim to reproduce the interstellar 
mass spectrum of nearby galaxies. We illustrate in Fig. 2 
the mass spectrum obtained with FF= 0.23 (red curve) in 
a shell geometry defined with i?min = 49 and i?max 64 
cells. The slope of the mass spectrum does not change no- 
ticeably decreasing the covering factor CF (ie. decreasing 

-^min ) ■ 

Let's mention that some models with other mass spectra 
(i.e. other filling factors FF) have been studied, such as 
p{in) ss m^^ ™ and p{m) ~ m^^'^'^. However, no notable 
change is found in any of our results changing only the mass 
spectrum in clumpy shell structures. 

2.2.3. Input spectra 

In the region close to Lya we assume that the spectrum con- 
sists of a flat UV continuum (i.e. constant in number of pho- 
tons per frequency interval) plus the Lya line, characterised 
by a Gaussian with an equivalent width EWint(Lya) and 
Full Width at Half Maximum FWHMi„t(Lya). All pho- 
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I 

7. 



FF = 0.07-p(m) = m*"' 
FF = 0.13 -p(m) = m"' 
FF = 0.23 - p(m) = m = " 




Number of cells in clumps 

Fig. 2. Variation of clump size adopting different filling factors 
FF in a shell geometry defined with -Rmin = 49 and -Rmax = 64 
cells. The filling factors used here are: FF= 0.07 (circles), FF = 
0.13 (diamonds) and FF = 0.23 (triangles). The mass spectrum 
obtained with FF =0.23 {p{m) cc m"^"'*) is the most consistent 
with observations of diffuse interstellar clouds. We thus adopt 
FF = 0.23 throughout this present study. 



tons are isotropically emitted from the center of our shell 
geometries. 

Throughout this work, we adopt FVFiJMi„f (Lya) — 
100 km s~^, as a typical value for the intrinsic width of the 
H recombination lines emitted in the ionized gas, observed 
in both nearby and distant starburst galaxies. Indeed, based 
on observations of the velocity dispersion in starburst galax- 
ies, this line width is comparable to the values measured 
from the velocity dispersion of CO and Ha lines in the star- 
burst galaxy cB58 (Tephtz et al. 2000; Baker et al. 2004), 
the velocity dispersion measured in several starbursts at z 
~ 2 by Erb et al. (2003) and in SMM J2135-0102 at z = 
2.32 by Swinbank et al. (2011). Furthermore, different val- 
ues of EWint(Lyct) , specified below if necessary, have been 
adopted. 

2.2.4. Output parameters 

For the present paper we are interested in the following 
quantities predicted by our Monte Carlo simulations: the 
average Lya escape fraction /osc(LyQ!), the Lya line pro- 
file and the escape fraction of continuum photons at UV 
and other wavelengths /esc (A). From these we also derive 
the observed Lya equivalent width i?Wobs(Lya), and the 
colour excess E{B — V). All quantities are computed from 
spatially integrating all the photons escaping our spheri- 
cally symmetric shells. They therefore correspond to aver- 
age properties for our homogeneous and clumpy structures. 

In the present study, the Lya escape fraction is com- 
puted from 



/osc(Lya) 



/!L/e^c(A)0(A)dA 



(5) 



where /esc (A) is the monochromatic escape fraction com- 
puted for typically 1000-2000 frequency points around the 
line center with a spacing of 20-10 km s~^, and 0(A) de- 
scribes the input line profile. The predicted Lya line profile 



can be computed a posteriori from our simulations for ar- 
bitrary input spectra (line -t- continuum), as described in 
Verhamme et al. (2006). /esc (Lya) is slightly dependent on 
the FWHM of the input line profile, but independent on 
the value of i?Wi„t (Lya). 

The UV continuum escape fraction /osc(UV) is com- 
puted redward of the Lya line. Assuming the Calzetti et al. 
(2000) attenuation law we can compute the corresponding 
colour excess as 



E{B-V) 



-2.5 



Calzetti 



fc(1235) 



log (/esc (UV)), 



(6) 



with fc(1235) = 11.4 according to the Calzetti law. 

From the escape fraction of radiation at the optical 
wavelengths listed in Table 2 we can also determine the 
true colour excess 



E{B - V), 



■cal 



-2.5 log 



fescjV) 
fesc{B) 



(7) 



In practice this is done by calculating the continuum escape 
fractions at 4350 A for the B band and at 5500 A for the 
V band (Table 2). 

The observed Lya equivalent width and the intrinsic one 
(i.e. input value of the source, before radiation transfer) are 
related by 



EWobs(Lya) /esc (Lya) 



(8) 



EWi„t(Lya) /esc(UV) ' 
where EWmti^^yct) is the intrinsic Lya equivalent width. 

2.2.5. Validation 

To test the radiation transfer, we have compared our re- 
sults to Witt & Gordon (2000), whose dust parameters are 
adopted in our calculations. We have constructed a clumpy 
shell model using the same discretisation and input parame- 
ters. The derived mass spectrum is in good agreement with 
these authors. The resulting continuum escape fractions, 
shown in Fig. 3, and other results are found in excellent 
agreement with Witt & Gordon (2000), which vahdates our 
code. 

3. The Lya and UV continuum radiation transfer 
in homogeneous and clumpy media 

In this section we study the radiative transfer of Lya and 
the UV continuum photons in homogeneous and clumpy 
shell geometries. 

3.1. The UV continuum escape fraction 

In Figs. 4 and 9, we examine the evolution of the UV 
escape fraction /esc(UV) in homogeneous and clumpy 
media. As shown in both figures, /esc(UV) depends on 
three main parameters: 

1) the dust content (tq) 

2) the dumpiness of the dust distribution, assumed to 
trace the Hi distribution (nic/ric) 

3) the covering factor (CF) 
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Fig. 3. Comparison of tiie UV-to-optical continuum escape frac- 
tion derived from our code with the results of Witt & Gordon 
(2000). The clumpy shell geometry studied here is built on a 
cartesian grid of A'' = 30'' cells and is characterized by the fol- 
lowing parameters: FF = 0.15, nic/nc = 0.01, i?min ~ 5 and 
^?max = 15. The evolution of the escape fraction at 25 different 
wavelengths is shown, from A = 0.1 /im to A = 30 nm, assum- 
ing three different dust optical depths tv in the shell structure 
(measured in V band): tv ~ 2, 4.5 and 8. The escape fractions 
obtained with MGLyo? are marked in red, whereas those ob- 
tained by Witt & Gordon (2000) are marked in grey. We find a 
very good agreement with Witt & Gordon (2000). 



The other parameters describing our shell geometries, 
the outflows (vexp), the Hi column density (iVni) and the 
temperature of the matter (b) do not show any effect on 
the UV escape fraction, as expected. 

Figure 4 illustrates the effects produced by both the 
dust content (i.e. t^) and the dumpiness of the dust 
distribution (i.e. njc/nc) on /csc(UV). We adopt here 
three media with the following conditions: FF = 0.23, CF 
= 0.997 and njc/nc = 1-0, 0.01, 0. Qualitatively, we can 
summarize the effects produced by both Tq" and njc/nc 
on /esc(UV) in the following way: 

- Ta'. in homogeneous and clumpy media, an increase 
in the dust optical depth Tq always produces a decrease in 
the UV escape fraction. 

- nic/n-c- a decrease in nic/nc from 1 to (i.e. from a 
homogeneous to an extremely clumpy dust distribution) 
always increases the UV escape fraction. A clumpy dust 
distribution produces indeed higher UV escape fractions 
compared with an equivalent homogeneous distribution of 
equal dust content (i.e. equal r^). 

Clumpy media are thus more transparent to UV con- 
tinuum radiation, as previously shown by Boisse (1990); 
Hobson & Scheuer (1993); Witt & Gordon (1996, 2000). 
The facts that the dust content concentrates in clumps 
and that the interclump medium becomes more optically 
thin, allow UV photons to escape any clumpy media in 
two different ways (Witt & Gordon 1996): first, like in 
homogeneous dusty media, UV photons have to scatter 
against few dust grains before escaping clumpy media (dust 
localized in clumps or in between clumps). But, in clumpy 
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Fig. 4. Evolution of the UV escape fraction /osc(UV) as a func- 
tion of the dust optical depth (r^) in three different shell geome- 
tries: a homogeneous shell {njc/nc = 1.00), a weakly clumpy 
shell {uic/nc = 0.01) and an extremely clumpy shell {nic/nc 
= 0.00). The clumpy media are defined with FF = 0.23 and CF 
— 0.997 (-Rmin ~ 49 and iimax ~ 64 cells). For comparison we 
also plot the curves for f^sc = exp(— tJ;) and /osc = exp(— r^), 
which describe respectively the upper and the lower limit of the 
continuum escape fraction in homogeneous media. 



media, UV photons take advantage of the weak opacity 
of the inter-clump medium, which allows them to escape 
more easily clumpy media than any homogeneous dusty 
geometry. Second, continuum photons can also directly 
escape clumpy media if several free spaces appear between 
clumps. However, that can be only possible in extremely 
clumpy media {njc/nc ~ 0). In this case, continuum 
photons are not affected by the weak dust content localized 
between clumps and can directly escape clumpy media 
getting through holes which appear between clumps. 

The covering factor CF is thus an important parameter 
controlling the UV escape fraction in clumpy media. Figure 
9 shows this dependence in the particular case of extremely 
clumpy shell geometries {njc/nc — 0). As expected, the 
UV escape fraction /osc(UV) decreases when the covering 
factor increases to unity (i.e. all lines-of-sight are covered by 
one or more clumps from the photon source when CF = 1). 
Furthermore, in the particular case of extremely clumpy 
shell geometries {njc/nc — 0), we can also notice that 
the covering factor CF provides a general lower limit for 
/esc(UV). As shown in Figs. 4 and 9, the UV escape frac- 
tion always converges on an asymptote /esc(UV) = 1 — CF, 
corresponding to the direct escape fraction. 

Besides this qualitative approach concerning the depen- 
dence of /osc(UV) to Ta, nic/nc and CF, Fig. 4 illustrates 
other quantitative results: in homogeneous geometries the 
UV escape fraction decreases very rapidly with the dust op- 
tical depth Tq. If we define /esc(UV) = e~'^'=^^ , the effective 
optical depth is equal to in the absence of scattering. 
With scattering the effective absorption increases, and one 
has Teff > Ta. In clumpy media the situation is different as 
photons can escape more easily, hence Teff < Tq. 

The escape fraction of the optical continuum photons 
evolves in the same way as for the UV photons in homo- 
geneous and clumpy media. Combining the escape fraction 
of both the B and the V-band in the same media as those 
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Fig. 5. Schematic representation of the way of propagation of 
UV continuum (red) and Lya photons (blue) in the "low con- 
trast" regime. The medium illustrated here is composed of high 
density clumps (containing Hi + dust), distributed in an inter- 
clump medium of low density (Hl+dust). In the "low contrast" 
regime, the Hi content in and in between clumps are relatively 
high, which renders both regions optically thick for the Lya pho- 
tons. In this way, the Lya photons can only escape the medium 
after undergoing multiple resonant scattering against Hi atoms, 
increasing their probability being absorbed by the dust. The 
UV photons are not affected by the presence of Hi atoms and 
propagate directly through the medium. 
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Fig. 6. Evolution of the Lya (blue lines) and UV continuum 
(red) escape fraction as a function of the dust optical depth (t^) 
for a homogeneous (triangles), a weakly clumpy (squares), and 
an extremely clumpy (spheres) medium. In these simulations, 
we adopt FF = 0.23, CF = 0.997 (iimin = 49 and ii^ax = 64 
cells), iVm = 10^^ cm"^ , Vexp = km s"^ and b= 40 km s~\ 
For comparison, the UV escape fractions (Fig. 4) are overlaid 
(red lines). For the extremely clumpy medium we note that for 
Ta larger than a certain limit Tc (here roughly equal to 3.5), 
/esc(UV) drops below /esc (Lya). 



studied in Fig. 4, we illustrate in Fig. 17 the evolution of 
the derived colour excess E(B-V). This figure can be used 
to translate the dust optical depth Tq of Fig. 4 in terms of 
colour excess E(B-V). 

3.2. The Lya radiative transfer in homogeneous and clumpy 
media: two regimes appear 

Besides the three main parameters which control the radia- 
tive transfer of the UV continuum photons (i.e. t^, njc/nc 
and CF), three other parameters also determine the radia- 
tive transfer of the resonant scattered Lya photons, namely 
Vexp, -^Hi and b. We now discuss the influence of these pa- 
rameters on the UV continuum and on Lya. For simplicity 
we here assume a constant value of b in all cells (clump 
or interclump). Overall we find that wc can identify two 
regimes where the Lya propagation is quantitatively dif- 
ferent, which we now explain. The separation between the 
regimes will be discussed after that (Sect. 3.3). 

3.2.1. The "low contrast" regime: homogeneous and weakly 
clumpy media 

Propagation of Lya photons in the "low contrast" 
regime 

We show in Fig. 5 the typical way Lya photons propagate 
in the "low contrast" regime. We deduce this propagation 
from our numerical simulations, studying the number and 
the location of each interaction between the Lya photons 
and the HI atoms in clumpy media. In this regime, we 
notice that the Lya radiative transfer is characterised by 
a (pseudo-)random walk in the medium, both in and in 
between clumps. 



Lya escape fraction /esc (Lya): 

Figure 6 shows the dependence of /osc(Lya) on the 
dust content (ja) of the medium, for different values of 
dumpiness {njc/nc)- The "low contrast" regime includes 
all curves of this figure, except the particular case njc/nc 
= 0.00 which belongs to the "high contrast" regime. Note 
that the quantity Tq can be related to the colour excess 
through Fig. 17. 

Qualitatively, we see that /esc (Lya) always decreases 
with increasing as expected. The same is true for increas- 
ing nic/nc- However, we see that for a given value of Ta, 
increasing njc/nc results in a faster decrease of /esc (Lya) 
than /esc(UV), the reason being the highly increased path 
length of Lya photons due to resonant scattering. Thus, in 
the "low contrast" regime Lya radiation is more vulnerable 
to dust than UV continuum radiation. 

A change of the covering factor CF has also a notice- 
able effect on /esc (Lya) in the "low contrast" regime (CF 
measuring the proportion of holes which appear between 
clumps). As the clumps cover an increasing fraction of the 
sky, it becomes indeed increasingly difficult for the photons 
to escape, and when CF w 1, /esc (Lya) drops drastically. 

Finally, the effect of a change of Vexp and iVni on 
/esc (Lya) is shown in Fig. 7. We notice that /esc (Lya) al- 
ways increases with increasing Vexp as well as with A'hi de- 
creasing. Since the effect of the expansion velocity is to shift 
the Lya photons away from the line center, they undergo 
progressively fewer scattering as Vcxp increases. In fact, for 
an intrinsic Lya line width of FW H Mint (Lya) (this value 
is 100 km s~^ in Fig. 7), a galactic outflow showing a ve- 
locity Vexp^ 2 X FWHMinti^'ya:) is enough to allow Lya 
photons escape as easil y as UV continuum photons. As ex- 
pected, an increase in A'hi always leads to a decrease in 
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Fig. 7. Evolution of the Lya escape fraction /esc (Lya) as a func- 
tion of dust optical depth Ta for various values of expa nsion 
velocity Vcxp and mean neutral hydrogen column density A^'hi- 
The clumpy medium studied here is built with: FF = 0.23, CF 
= 0.997 (i?min = 49 and i?max = 64 cells), ruc/ric = 0.01 and 
b= 40 km s~^. The expansion velocities Vexp and Hi column 
densities that we adopt are: Vexp = 0, 100 and 250 km s~^ and 
iVm = 10^^ cm"^ and 2x10^° cm"^ . It is seen that /osc(LyQ) 
increases with increasing Vcxp,as well as with decreasing A^hi- 

/csc(Ly«), since more neutral hydrogen implies more scat- 
terings, and hence an increased total path length before 
escape, resulting in an increased probability of being ab- 
sorbed. 

Quantitatively, both Figs. 6 and 7 allow us to general- 
ize the fact that, in the "low contrast" regime, we always 
obtain: 

/esc (Lya) < /osc(UV) (9) 

The strict equality /esc (Lya) = /esc(UV) is only met when 
Lya photons are able to avoid scatterings altogether, i.e. 
for a sufficiently high expansion velocity, or for a very low 
Hi column density. 

Lya equivalent width EW(Lya): 

Combining the definition of the Lya equivalent width (Eq. 
8) and Eq. 9, we always have: 

EWobsiLya) < EW,nt{Lya) (10) 

In the "low contrast" regime, the observed Lya equivalent 
width EWohs{Lya) is thus always lower or equal to the in- 
trinsic one EWintiLya). In other words the Lya equivalent 
width is not "boosted" by clumping in this regime. 

3.2.2. The "high contrast" regime: extremely clumpy shell 
geometries 

The "high contrast" regime of the Lya radiative transfer is 
defined by clumpy media showing a very low ratio njc /tiq , 
i.e. a high density contrast at least njc/nc ^ 1.5 x 10^^ 
for the input parameters considered throughout this 
study (see Sect. 3.3). To describe the main effects and 
peculiarities of this regime we here restrict ourselves to the 
most extreme case with nic/nc = 0. 




Higli deosil^' : HI + dust iDter-cliimp medium : emptj' 



Fig. 8. Schematic representation of the way of propagation of 
the UV continuum (red) and the Lya (blue) photons in the "high 
contrast" regime. The medium illustrated here is composed by 
high density clumps (Hi -f dust), distributed in an empty in- 
terclump medium. A fraction of the Lya photons scatter on the 
Hi atom comprising the surface of the clumps. The rest of the 
Lya photons, as well as the UV continuum photons, pierce the 
clumps where they may be absorbed by dust. 



Propagation of Lya photons in the "high con- 
trast" regime: 

The way Lya photons propagate now differs qualita- 
tively from the radiative transfer in the "low contrast" 
regime. The details of the way we deduce the propagation 
of Lya photons in the "high contrast" regime is given in 
Sect. 4, where we study the Lya line shape. We sketch in 
Fig. 8 the propagation of Lya photons in the "high con- 
trast" regime. The Hi content distributed between clumps 
is now weak enough that scattering between clumps can 
be neglected, and for a fraction of the Lya photons, the 
radiative transfer is characterised by rebounces on the 
clumps, as originally suggested by Neufeld (1991). The 
remaining Lya photons propagate in the same way as UV 
photons, that is penetrating the clumps, being exposed to 
the dust, or escape the medium freely (if CF < 1). 

Lya escape fraction /esc (Lya): 

Figure 6 shows the evolution of /esc (Lya) as a func- 
tion of the dust content (i.e. t^) in the "high contrast" 
regime (jiic/nc = 0). Again, note that the quantity fa 
can be related to the colour excess through Fig. 17. From 
Fig. 6 it is evident that /esc (Lya) always decreases as 
increases. However, the decrease is slower than that of 
/esc(UV), allowing the curve of /esc (Lya) to cross that of 
/osc(UV) at a certain optical depth Tc (tc « 3.8 in Fig. 
6). This is not possible in the "low contrast" regime, and 
it is this quantitative difference that defines the threshold 
between the low and the high contrast regime. 

Figure 9 illustrates the effect of a change of the cover- 
ing factor CF on /esc (Lya) in the "high contrast" regime. 
Again, we see that both /esc (Lya) and /csc(UV) always 
decreases with increasing CF. However, whereas /osc(UV) 
approaches asymptotically the value 1— CF (corresponding 
to all clumps being fully opaque to the UV so that escape is 
possible only through direct escape), /esc (Lya) maintains 
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Fig. 9. Lya (blue lines) and UV (red lines) escape fraction as a 
function the dust optical depth tJ" for different values of covering 
factor CF in the "high contrast" regime {njc/nc = 0.00). In this 
figure, all media have for common parameters: FF = 0.23, nic/ric 
= 0.00, iVm = 10^^ cm"^ Vexp = km s"^ and & = 40 km s"\ 
The values Tci, Tc2 and Tcs, corresponding to the values of tJ" 
where the Lya and UV escape fractions of CF\, CF2 and CF^ 
cross, are marked. While /es5c(UV) always converges towards the 
limit /esc(UV) = 1 - CF in extremely clumpy media, /osc(LyQ!) 
shows higher values. 



a higher value. Furthermore, we can notice from the Fig. 
9 that the value of the critical optical depth Tc (where 
/escCLya) crosses /csc(UV)) strongly decreases as CF de- 
creases. 

The effect of a change of Voxp and A^hi on /csclLyct) is 
shown in Fig. 10. We can see that these effects are different 
and more complex than those of the "low contrast" regime. 
More precisely, two different domains appear in the "high 
contrast" regime, below and above the critical optical depth 
Tc (where /osc(Lya) crosses /csc(UV)). Below Tc, we notice 
that /esc(Lya) increases with increasing Vcxp as well as 
with A^Hi decreasing. Such behavior is the same than those 
of the "low contrast" regime. But above Tc, the opposite 
effect is observed, where an increase of /csclLya) results 
from a decrease of Vcxp and an increase of A^hi- These 
distinct behaviors can be understood as follows (see Fig. 8): 

- Ta ^ Tci In this domain, Lya photons are more 
vulnerable to dust than UV photons. The dust content 
being relatively low, each clump is optically thin for UV 
radiation which allows UV photons to escape directly 
the medium getting through the clumps. However, Lya 
photons have to scatter off of the surface of a high number 
of clumps before escaping (Fig. 8), which increases the 
probability of being absorbed by the dust. An increase of 
Voxp increases the Lya escape fraction /osc(LyQ;). Indeed, 
increasing the expansion velocity Vcxp, all Lya photons are 
Doppler shifted out of resonance, forcing them to pierce 
the clumps and thus to escape the medium a s eas ily as UV 
photons. For the same reason, a decrease of Nhi increases 
/esc (Lya) because decreasing the Hi density in clumps. It 
renders therefore clumps more transparent to Lya photons. 

- Ta Tc'. In this domain, Lya photons are less af- 
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Fig. 10. LyQ (blue lines) and UV (red line) escape fractions as 
a function of the dust optical depth Ta in the "high contrast" 
regime, for various values of expansion velocity Voxp and hy- 
drogen column density A^'hi. The extreme clumpy media have 
in common: FF = 0.23, CF = 0.997 (i?min = 49 and ii^ax = 
64 cells), nic/ric — 0.00 and 6 = 40 km s"'^. We adopt as well 
the following expansion velocities Vexp and Hi column densities 
iVm: Vexp = 0, 100 and 250 km s"^ and iV^^^ = lO^'' cm"^ and 
2 X 10^° cm~^. We indicate the critical dust optical depth Tc 
where the curve of /eac(Lya) for Vexp = km s~^ cross the 
curve of /csc(UV). 

fected by the dust than UV photons. Indeed, while UV 
photons are now strongly absorbed by the high dust 
content embedded in clumps, Lya photons can avoid inter- 
action with dust scattering off of the surfaces of clumps. 
An increase of Voxp results in a decrease in /esc (Lya). 
Increasing the expansion velocity all Lya photons are 
Doppler shifted out of resonance, forcing them to pierce 
the clumps, where they are strongly absorbed by the dust 
like UV photons. For the same reason, /esc (Lya) now 
increases with increasing TVhi because increasing the Hi 
content of the clumps the probability of Lya photons 
scattering off the clumps increases, having the journey 
confined to the dust-free interclump medium. 

Finally, as in the "low contrast" regime, we can no- 
tice that a velocity Vgxp ^ 2 x FW H MintO^yct) is enough 
to prevent any scattering on the clumps. In this case, Lya 
photons escape the medium in the same way than UV 
photons. 

Lya equivalent v^^idth EW(Lya): 

Qualitatively, Figs. 6, 9, and 10 reveal that in the 
"high contrast" regime /esc (Lya) can be both higher or 
lower than /esc(UV), depending on the actual value of Tq. 
From the definition of Tc (i.e. the value of t^ where the 
curves of /esc (Lya) and /esc(UV) cross each other), we 
have 

/csc(Lya) < /esc(UV) if < Tc, (11) 
and 

/csc(Lya) > /esc(UV) if t^ > Tc, (12) 
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where we note that Tc is mainly a function of CF. 

Combining the definition of the Lya equivalent width 
(Eq. 8) with both Eqs. 11 and 12, we obtain: 

EWobs{Lya) < EW,nt{Lya) for Ta < (13) 

and 

EWobs{Lya) > EW„,t{Lya) for > t,, (14) 

The fact that /osc(Lya) can exceed /osc(UV) thus allows 
for an enhancement ("boost") of the equivalent width of 
the Lya line. However, as summarized in Sect. 5.2, such an 
enhancement can only occur under strict physical condi- 
tions, concerning the kinematics (i.e. an expansion velocity 
Vexp ^ 2 X FW H AIint(Lya)) , the dumpiness and the dust 
content of the clumpy ISM. 

3.3. The critical ratio nic/ric separating the two regimes of 
the Lya radiative transfer 

We can quantify the distinction between the two regimes of 
the Lya radiative transfer by the low to high density ratio 
njc/nc. Figure 11 illustrates this limit njc/nc as a func- 
tion of the average Hi column density TVhi in our clumpy 
media. This limit is defined in the following way. Above the 
curves shown in Fig. 11, it is impossible to obtain /osc(Lya) 
> /esc(UV) in models with physical parameters listed in 
Table 3. Such media belong to the "low contrast" regime. 
Conversely, the area localized below the curves corresponds 
to the "high contrast" regime where it is possible to observe 
the inequahty /esc (Lya) > /osc(UV). 

All points shown in Fig. 11 have been obtained study- 
ing a clumpy medium defined with: FF = 0.23, CF — 0.997, 
Vexp = km s~^, Tq = 25 and 6 = 40 km s^^. All media built 
with different parameters show a critical ratio nic/nc (at 
the limit between the low and the high contrast regimes) 
lower than that shown in Fig. 11. Each curves shown in 
this figure corresponds to three different physical conditions 
applied in our clumpy media. The star dots are obtained 
assuming an unique turbulent velocity 6 = 40 km s^^ in 
and in between clumps, and an inter-clump medium com- 
posed by both Hi atoms and dust grains. The circle dots 
are obtained assuming the same turbulent velocity (5) but 
a dust-free interclump medium. Finally, the squares are ob- 
tained applying two different temperatures in clumps and 
in between clumps {b = 40 km s~^ in clumps and T = lO^K 
between clumps) and assuming an interclump medium only 
composed by Hi atoms. This case is discussed in Sect. 3.4. 

From the single temperature case (star and circle dots) 
we can already mention four main results concerning the 
border separating the two regimes of the Lya radiative 
transfer: 

The critical ratio njc/nc separating both regimes 
is very low: studying a large range of Hi column density 
■^Hi [10^'', 10^^] cm~^, we notice that the limit separating 
both regimes is reached for very weak ratios njc/nc 
([1.5 X lO""*, 1.3 X 10"^]). That suggests that the "high 
contrast" regime can only be found in galaxies showing 
the most extremely clumpy ISMs (composed only by cold 
clouds of neutral hydrogen gas embedded in an extremely 
ionized interclump medium). The critical ratio has to 
decrease with increasing A'hi to maintain a sufficiently low 
column density between the clumps. 



0|^_^_ ^ — . . 

J L )l\ One temperature - inteiclump : HI + dust 

One temperature - intercltmip : HI 
-2 - ■ T^wo temperatures - interclump : HI 

-3 r 

■ w « low contrast » regime 
-4 r * « 

- ■ X 

■ K 



« high contrast » regime 
-8 - - 
.9 r ■ I , , ,_] , . J . , J , , J , , .J ._ 

17 18 19 20 21 22 

log(NHIJ (cm -) 

Fig. 11. Evolution of the critical ratio nic/nc separating the 
"low contrast" regime to the "high contrast" regime of the Lya 
radiative transfer in clumpy shell geometries. For highest ratio 
"■/c/wc, all clumpy media belong to the "low contrast" regime, 
whereas for lower nic/ric the "high contrast" regime is ob- 
served. Three curves are represented in this Fig. corresponding 
to three different physical conditions applied to our clumpy me- 
dia. The star and circle dots have been obtained assuming 1) a 
single temperature (6 = 40 km s"'^) in and in between clumps, 
but 2) not the same composition between clumps (Hi -I- dust for 
the star dots and only Hi atoms for the circle dots). The square 
dots assume rather 1) a warmer temperature between clumps 
than in clumps (T = 10*' K between clumps and 6 = 40 km s~^ 
in clumps) and 2) an inter-clump medium composed only of Hi 
atoms. 



The interclump medium can be optically thick 
for Lya photons on the border separating the two 
regimes: In Fig. 12 we show the limit of Fig. 11, but 
translated in terms of HI column density between clumps 
-^Hi.ic = (1 ~ FF) xn/cL. While the interclump medium 
is optically thick for Lya photons at line center above 
A^Hi.ic — 6 xlO^^ cm~^ (assuming 6 = 40 km s~^ between 
clumps)^, we can clearly see that the "high contrast" 
regime can extend to somewhat higher interclump column 
densities, into the optically thick regime. Nevertheless, 
this is only possible if the temperatures in and in between 
clumps are the same (star and circle dots). If the temper- 
atures in and in between the clumps differ (square dots), 
the interclump medium has to be optically thin in order to 
observe the "high contrast" regime (see Sect. 3.4). 

The limit njc/nc separating both regimes is 
not affected by the presence of dust between 
clumps: The critical ratio nic/nc represented by both the 
star and the circle dots on Fig. 11 are the same. Thus, the 
presence of dust between clumps has no effect on the limit 
separating both regimes of the Lya radiative transfer. 



^ The optical depth at line center is ro — 3.31 x 
^Q-i4jv-i/2^^^ = 3.31 X 10""(12.85kms-V^)^m> e.g. 
Verhamme et al. (2006). 
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Fig. 12. Same as in Fig. 11, but showing the evolution of the Hi 
column density between clumps A*'hi,ic (= (l-FF)xn/cL, with 
L as the physical size of the shell), as a function of the total 
A'hi- When the temperature in and in between clumps are the 
same (star and the circle dots), the "high contrast" regime can 
be observed applying an optically thick interclump medium for 
Lya photons (ie. from A^hi.ic > 6x10^'' cm^^ with fe = 40 km 
s~^). But, assuming a higher temperature between clumps than 
those of the clumps (square dots), the "high contrast" regime 
disappears as the interclump medium becomes optically thick for 
Lya photons (ie. from iVmjc > 3x10^* cm~^ for T = lO'^K). 

3.4. The effects of an in homogeneous temperature on the 
Lya radiative transfer 

If we want to render the physics of our clumpy media 
more realistic, we must assume different temperatures in- 
side and between the clumps. There is ample evidence for 
such a multi-phase ISM. For example, for the stability of 
the clumps a pressure equilibrium should intervene between 
the two phases of our clumpy media, implying a higher tem- 
perature in the interclump medium. Also, it is well known 
that different regions coexist in the real ISM of any galaxies 
(McKee & Ostriker 1977). In particular, we can distinguish 
the warm neutral atomic medium (WNM, where Hi atoms 
are present in the atomic form) to the hot ionized medium 
(HIM, where Hi atoms are in great majority ionized). To 
explore in a simplified manner these effects we have made 
calculations assuming a temperature of T = 10® K in the 
interclump medium (as those measured in the HIM), but 
a lower temperature in all clumps {b = 40 km/s, by anal- 
ogy to the WNM) . We then examine the effects on the Lya 
radiative transfer and on the ratio njc/nc separating the 
two regimes identified above. 

Qualitatively, the Lya radiative transfer properties be- 
have in a similar fashion for different interclump tempera- 
tures. However, as shown in Fig. 11, the ratio njc/nc sepa- 
rating both regimes is found to be lower than for the case of 
constant temperature. In other words, the "high contrast" 
regime is more limited when the temperature in the inter- 
clumps medium is higher than those in clumps. The reason 
for this is the following: on one hand the optical depth in 
the interclump medium decreases (with T~^/^), simulating 
thus a medium with even lower density, i.e. a higher con- 
trast. On the other hand, the temperature increase leads 
to a larger frequency redistribution of the scattered Lya 
photons, which eases their escape due to higher frequency 
shifts. This effect dominates over the former, rendering thus 



the clumps more transparent to Lya radiation, where they 
are strongly absorbed by the dust. This explains why even 
higher density contrasts are needed to achieve significant 
Lya "rebounce" on the clumps, if the interclump medium 
is hotter than the clumps. 

In terms of interclump column densities the limit be- 
tween the regimes is shown in Fig. 12. In contrast to the 
case of uniform "cold" temperatures, the limit is now found 
at quite low column densities of the interclump medium, 
corresponding to an optically thin regime. Indeed, such low 
column densities are needed if one wants to avoid significant 
scattering of Lya with the corresponding high frequency 
shifts in a hot interclump medium. 

4. Lya line profiles formation in homogeneous and 
clumpy shells 

In this section we give an overview of the different emergent 
Lya line profiles produced in the expanding homogeneous 
and clumpy shell geometries of our model. 

4.1. Lya line profiles from dust free homogeneous and 
clumpy shell geometries 

First we shall consider the case of a dust free ISM and 
examine how the Lya profiles are modified by a clumpy 
ISM structure. The line profiles shown in this section are 
obtained assuming an intrinsic Lya line characterized by 
EW„,t{l^ya) = 80 A and FWHMint = 100 km s'^. 

We know that in dust free cases the total Lya flux is 
preserved, and since the continuum is not attenuated (due 
to the absence of absorption), the Lya equivalent width is 
thus preserved, in other words the observed EW is identical 
to the intrinsic one. This holds obviously both for homo- 
geneous and clumpy structures. The only effect of clumps 
is to modify the exact frequency redistribution of photons, 
i.e. the shape of the emergent Lya line profile. However, as 
we will see, the changes to the line profile are only relatively 
small, when the covering factor of the clumps is large. 

We first examine the Lya line profiles for dust free ho- 
mogeneous and clumpy structures with low density con- 
trast (i.e. clumpy structures showing an optically thick in- 
terclump medium for Lya photons). In Fig. 13 we study 
such structures built with the following parameters: iVni — 
2 X 10^0 cm-2 , & = 40 km s'^ and Vexp = 0, 100, 300, 400 
km s^^. It clearly appears that the Lya line profiles emerg- 
ing from dust free weakly clumpy shell geometries do not 
show any noticeable difference compared to homogeneous 
structures with the same/corresponding properties. 

In homogeneous and weakly clumpy media, the mecha- 
nisms of formation of the line profiles, as well as their de- 
pendance on the parameters Vexp, -^hi and b, are identical 
to those explained in detail in Verhamme et al. (2006) for 
homogeneous shell structures. In other words, in the dust 
free case, weakly clumpy media do not significantly differ 
from the homogeneous ones in terms of line profiles. 

Turning now to clumpy media with large density con- 
trasts (i.e. clumpy structures showing an optically thin in- 
terclump medium for Lya photons) we find again very 
similar line profiles as in the homogeneous case, as also 
shown in Fig. 13. Compared with the line profiles observed 
from static weakly clumpy and homogeneous media we now 
see a central peak at line center (uobs~0 km/s). The for- 
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Fig. 13. Comparison between the Lya line profiles emerging 
dust free homogeneous shell (dotted blue lines), weakly clumpy 
shell (dashed magenta lines) and extremely clumpy shells (thick 
red lines) under different expansion velocities Vcxp- All homoge- 
neoues and clumpy shells have for common parameters: A''hi = 
2x10^° /cm-2 , b= 40 km s"\ ^ = and Vexp = 0, 100, 300, 
400 km s~^. Furthermore, both weakly and extremely clumpy 
media are built adopting the following parameters: FF = 0.23, 
CF = 0.997 {Rmin = 49 and Rmax = 64 cells) and iiic/nc = 
0.01 (weakly clumpy shell) and 0.00 (extremely clumpy shell). 



mation of this central peak is indeed made possible when 
the interclump medium is optically thin for Lya photons. 
These photons can then propagate in two different ways 
in the interclump medium: either by escaping through the 
holes which appear between clumps, or scattering off of the 
surface of clumps. Both features preserve the intrinsic fre- 
quency of the Lya photons, which allows to form the central 
peak seen at center of the Lya line (wobs~0). The covering 
factor governs mostly the importance of the central emis- 
sion, as shown in Fig. 14 for a static shell. As expected, 
the central emission increases with decreasing the covering 
factor. 

Besides the predicted Lya line profiles in dust-free 
clumpy shell geometries behave in the same way as already 
shown for homogeneous shells by Verhamme et al. (2006), 
with the main parameters determining the Lya profile be- 
ing the expansion velocity Voxp, the mean HI column den- 
sity iVui, and the Doppler parameter b. 
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Fig. 14. Effect of the covering factor CF on the Lya line pro- 
files in extremely clumpy shell geometries. We study here a static 
extremely clumpy shell geometry defined with the following pa- 
rameters: FF = 0.23, nic/nc = 0.00, iVm = 2x10^" /cm"^, 
Vexp = km s~^. To = and 6 = 40 km s~^. The only differ- 
ence between both panels concerns the covering factor CF of the 
clumpy shell geometries: CF = 0.99 (left panel) and CF — 0.95 
(right panel). We can notice that decreasing CF, the relative 
intensity of the central peak in the Lya line profile is increased. 
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Fig. 15. Effects of the dust optical depth Ta on the Lya line 
profiles emerging three different static shell geometries: a homo- 
geneous shell {nic/nc = 1.0 - top line), a weakly clumpy shell 
{nic/nc = 0.01 - middle line) and an extremely clumpy shell 
{nic/nc ~ 0.00 - bottom line). Both the weakly and the ex- 
tremely clumpy shell geometries are those studied, respectively, 
in figures 14 and 17. We thus adopt here the same physical con- 
ditions: FF = 0.23, CF = 0.997, Nm = 2x 10^° /cm'^ , v^xp 
= km s"'^ and & = 40 km s"'^. Note the variable scales of the 
different sub panels. An enhancement of EW(Lya) (produced 
by the Neufeld effect) occurs only for the physical conditions 
corresponding to the bottom right panel (tJ" = 3). 



4.2. Lya line profiles from dusty homogeneous and clumpy 
shell geometries 

We now examine the main effects produced by dust on the 
Lya line profiles emerging from clumpy shell geometries. 

In Fig. 15 we illustrate the evolution of the line profiles 
predicted for static homogeneous and clumpy shell geome- 
tries as a function of the dust optical depth Tq. The top line 
shows the homogeneous case, the middle the low density 
contrast, and the bottom line the clumpy medium with a 
high density contrast. Let us first examine the homogeneous 
and low density contrast cases {nic/nc = 1.00 and 0.01 re- 
spectively). Increasing from to 1, the Lya line still ap- 
pears in emission in both cases. But, a clear decrease in both 
the width and the intensity of each peak are noticed. For 



Tq = 3, more than 99.8 % and 98 % of the Lya photons are 
absorbed by the dust, respectively in the homogeneous and 
the clumpy media. Therefore, an absorption profile emerges 
from the homogeneous medium, while faint emission line is 
predicted from the weakly clumpy medium. Higher dust op- 
tical depths Ta are needed to obtain absorption line profiles 
from weakly clumpy media, typically fa ^ 35 in the case 
of nic/nc = 0.01. 

For extremely clumpy shell geometries (jiic/nc = 0.00, 
Fig. 15), the evolution of the Lya line profile is different. 
Increasing Ta from to 1, we notice a clear decrease of the 
width of both lateral peaks, as well as an increase of the rel- 
ative intensity of the central peak. The photons composing 
the central peak indeed interact very weakly with the dust. 
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Fig. 16. Same as Fig. 15 but for an expanding shell with Vcxp= 
300 km s~^. Note that no boost of Lya with respect to the 
continuum is found in any of these models, since the expansion 
velocity is too high. Discussion in the text. 

which explains why this peak becomes the dominant one in 
the line profile as increases. When is further increased 
from 1 to 3, both lateral peaks are destroyed by dust. The 
central peak thus becomes the only peak composing the 
line profile above ^ 3. It is interesting to note that we 
cannot obtain absorption line profiles in extremely clumpy 
media {njc/nc = 0.00). Indeed, as the Lya photons com- 
posing the central peak interact very weakly with dust, they 
are always able to escape the medium for any dust optical 
depth Ta, giving thus rise to an emission line. Comparing 
the relative escape of the Lya and UV continuum photons, 
we note that the only case in Fig. 15 where the Lya equiv- 
alent width is (slightly) enhanced is found in the bottom 
right panel. Indeed, in this case is close to critical dust 
optical depth Tc ~ 3 for this example of extremely clumpy 
medium, where we expect such an enhancement (cf. Sect. 
3.2.2). 

We now turn to a case with outfiows in Fig. 16. In this 
figure, we adopt otherwise identical parameters to those 
shown in Fig. 15. Increasing in any media (homoge- 
neous or clumpy), we first notice a quick decrease of the 
intensity of both the dominant red peak (those shifted at 
Vobs = 2xvexp) and the small blue bump. This is due to 
higher number of scatterings these photons undergo, which 
increases their destruction probability, as already discussed 
by Verhamme et al. (2006). For the highest dust content 
(tq = 3), the Lya line escaping homogeneous and clumpy 
media exhibits an asymmetric profile, but whose the dom- 
inant peak is found at line center (wobs=0). Finally, like in 
the static case, we notice that the intensity of the line in- 
creases as the dumpiness of the medium increases. In none 
of the cases shown here we find a " boost" of the Lya equiv- 
alent width, since the velocity is too large. 

5. Discussion 

5.1. Effects of a clumpy ISM on the radiation attenuation 

Given the evolution of the Lya and continuum escape frac- 
tion in homogeneous and clumpy systems (Sect. 3), it is 
clear that a clumpy medium always produces higher Lya 



and continuum escape fraction compared with an equiv- 
alent homogeneous medium of equal dust and hydrogen 
mass. This main result was demonstrated by several previ- 
ous studies focused on the transfer of the continuum radia- 
tion in clumpy media (Boisse 1990; Hobson & Scheuer 1993; 
Witt & Gordon 1996, 2000; Varosi & Dwek 1999), but also 
from other studies focused on the Lya line (Neufeld 1991; 
Hansen & Oh 2006). 

The attenuation of the radiation in a galaxy is thus 
strongly dependent on both the dust content and the dust 
distribution around the radiation sources. For illustration, 
we show in Fig. 17 the dependence of the colour excess E(B- 
V) on both the dust content (tq) and the dumpiness of the 
dust distribution in the shell geometries studied throughout 
this paper. In this figure we compare two different defini- 
tions of the colour excess: E{B — V)reah which corresponds 
to the exact colour excess because estimated from the origi- 
nal definition of the colour excess (from the V and B bands), 
and E{B — V)caizetti which is estimated from both the 
Calzetti attenuation law (Calzetti et al. 2000) and the UV 
escape fraction (see Eqs. 6 and 7). In practice, the Calzetti 
attenuation law is usually used to estimate the dust attenu- 
ation in starburst galaxies. It is then E{B—V)caizetti which 
would be measured by an observer. The Fig. 17 allows then 
us to see in which extend the colour excess E{B—V)caizeui, 
from the Calzetti law, deviates from the real colour excess 
E{B — V)reai as a function of Tq and the dumpiness of the 
dust distribution. 

In a general way, we can notice that the dumpiness 
of the dust distribution strongly affects both the colour 
excess E{B — V)reai and E{B — V)caizeti (Witt & Gordon 
2000). The colour excess decreases as the dust distribution 
is clumpy and as the dust optical depth decreases in media. 
Comparing now both definitions of E{B — V)caizetti and 
E{B — V)reai, we Can notice that E{B — V)caizetu does not 
reproduce very well the real evolution of the colour excess 
E{B — V)reai- This deviation between both definitions is 
explained by a clear evolution of the attenuation law (which 
measures, at each wavelength, the reduction in the stellar 
flux from a dusty ISM) as the dust distribution and the 
dust content change in media. As mentioned in Witt & 
Gordon (2000), this divergence shows that the use of the 
same and unique attenuation law in the analysis of a large 
sample of galaxies (which show different dust geometries 
and dust content) can become a source of error in the dust 
attenuation correction for individual galaxies. 

5.2. High Lya EWs and the Neufeld model 

5.2.1. Physical conditions needed in the ISM 

According to our study of the Lya transfer in clumpy 
media, there exists a regime in which the Neufeld scenario 
works. This regime corresponds to the "high contrast" 
regime, as explained in detail in Sect. 3.2.2. It is only 
found in the most extremely clumpy shell geometries of 
our model, that is composed by clouds of Hi and dust 
embedded in an interclump medium close to be optically 
thin for Lya photons. However, even in this conflguration, 
the Neufeld model only works when the following five main 
conditions concerning the dumpiness, the kinematic, the 
dust content and the spatial distribution of the clumps 
around the stars are fulfilled: 
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Fig. IT. Evolution of the colour excess E(B-V) as a function of 
Ta in the homogeneous and clumpy shell geometries studied in 
both sections 3 and 4. The shell geometries are defined by the 
following parameters: FF = 0.23 and CF = 0.997 (i?„i„ = 49 
and Rmax = 64 cells), nic/nc = 1.0, 0.01 and 0. Two differ- 
ent colour excess are shown in this figure. Firstly E{B — V)reai, 
which corresponds to the exact value of the colour excess be- 
cause estimated from the original definition of the colour excess 
(from the V and B bands). Secondly E{B — V)caizetu which is 
estimated from both the attenuation law of "Calzetti" (Calzetti 
1997) and the UV escape fraction. An observer would rather 
measure E{B — V)caizetti in practice. 



The galaxy outflow has to be relatively slow: 

Assuming an intrinsic Lya line width FW H Mint(LysL) , 
a galactic outflow with an expansion velocity Vexp ^ 2 x 
i^VFifMint (Lya) km s~^ is needed to be able to enhance 
EW(Lya) under the Neufeld scenario, hi starburst galax- 
ies, the width of the intrinsic Lya line is lower than 100 km 
s-i (TepUtz et al. 2000; Baker et al. 2004; Erb et al. 2003; 
McLinden et al. 2011), which implies an expansion velocity 
Voxp lower than 200 km s~^ in the ISM. We illustrate this 
limit in Fig. 18. This figure shows the evolution of the 
ratio /osc(Lya)//osc(UV) as a function of the dust content 
(measured here in terms of colours excess E(B — V)caizetti) 
in an extremely clumpy shell geometries. We here adopt 
b = 12.8 km s-\ TVm = 10^^ and 2 x 10^" cm'^ and 
Vexp = 0, 100, 200 km s^^, typical of values obtained in 
the analysis of high-z Lya line profiles (Verhamme et al. 
2008) . Finally, let us mention that the curves shown in this 
figure illustrate the highest enhancements of EW(Lya) 
we can obtain adopting such physical conditions in a 
clumpy medium; they reach up to a factor 3-4. Adopting 
FWHM{Lya)int = 100 km s^^ in Fig. 18 we can notice 
that no significant enhancement of EW(LyQ!) is obtained 
for Vcxp ^ 200 km s^^. Above such expansion velocity, 
all Lya photons are Doppler shifted out of resonance, 
preventing them to scatter off of the surface of clumps and 
to escape clumpy ISMs more easily than UV continuum 
photons. 

The galaxy outflow has to be relatively uni- 
form (constant velocity): Lya photons can scatter on 
the surface of clumps (Fig. 8) under the condition that 
each clumps move weakly each other. Should the opposite 



occur (that is assuming a random component Wrandom 
in the velocity of each clumps) a strong Doppler shift 
can occur between clumps which prevents Lya photons 
to scatter against clumps anymore. Finally, such effect 
strongly decreases the Lya escape fraction /osc(Lya) in a 
clumpy ISM. In Fig. 19, we illustrate how the enhance- 
ments of EW(Lya) shown in Fig. 18 are affected by a 
nonuniform outflow. In this figure, we assume a radial 
and random velocity Vdump for each clump, such that 

Wclump = Voxp + ^random, whcrC Wrandom = ^Wmax with 

r G [^1,1] a random number. We notice a clear decrease 
of the ratio /osc(Lya)//osc(UV) as ^random increases. 
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Fig. 18. Evolution of the ratio /esc(Lya)//esc(UV) (i.e. 
EWohB{Lya) / EWint[Lya)) as a function of the colour excess 
E{B — V)caizetti iu four extremely clumpy media. We apply some 
realistic physical conditions of LAEs to each media: nic jric ~ 

0. 00, FF = 0.23, CF = 0.997 (7?mi„ = 49 and Rmax = 64 cells), 
b = 12.85 km s"^ (i.e. T = 10^ K in clumps), IVm = 10^^ 2x 10^" 
cm~^ and Vexp = 0, 100, 200 km s~^. /csc(Lya) is derived as- 
suming an intrinsic line width FW H Mint = 100 km s~^. From 
the left to the right, the models are calculated for — 0, 0.3, 

1, 2, 5, 10, 15 and 25. 

The Hi content between clumps must be 
extremely small: The Neufeld model was originally 
developed assuming an interclump region sufficiently poor 
in Hi atoms, such that it is completely transparent to 
Lya photons. In our simulations, we have identified the 
allowed Hi content in the interclump medium to allow the 
Neufeld scenario to work. There is indeed a certain Hi 
limit above which Lya photons cannot freely propagate 
between clumps, preventing them to escape the medium 
more easily than UV photons. In all physical conditions, 
our simulations confirm that the Neufeld model only works 
if the interclump medium stays optically thin for Lya 
photons, that is if the radial Hi column density of the 
interclump medium {Nm ic) is lower than 3 x lO^'' cm~^ 
(with a temperature of T = 10^ K between clumps). For 
instance, focussing on the clumpy shell geometries studied 
in Fig. 18, we notice that no enhancement of EW(Lya) is 
obtained if the radial Hi column density between clumps 
exceeds 1.5 x lO^** cm^^ (for the curves A'hi = 10^^ 
cm-^) and 2.3 x 10" cm^^ (foj. ^j^g ^^^^^ = 2 x 10^° 
cm^^). In terms of ratio nic/nc, such limits correspond to 



14 



F. Duval , D. Schaerer , G. Ostlin , P. Laursen: Lya radiative transfer in a clumpy ISM 




0.1 D.2 0.3 0.4 0.5 



Fig. 19. Effect of a nonuniform outflow on the ratio 
/csc(LyQ)//Gsc(UV) in the same clumpy shell geometry stud- 
ied in Fig. 18 (FF = 0.23, CF = 0.997, iVm = 10^^ cm'^ and 
b = 12.8 km s"^). In this fi gure we assume a radial and random 
velocity ?;random for each clump with a maximum velocity of 
km s^^ (circle), 50 km s~^ (square) and 100 km s^^ (triangle). 



a density ratio of 6.90 x 10 ^ (for the curves A^hi = 10^^ 
cm-2) and 3.45 x 10^"^ (for the curve iVm = 2 x lO^" 
cm"^). In reality, lower densities ratios can be observed in 
a real ISM, if the cold clouds of neutral Hi (T = IC* K, 
n-m = 0.3 cm""^) are embedded in a very hot and ionized 
interclumps medium (T = 10^ K, nn ~ SxlO^'^ cm^^ and 
xh/ < lO^^ '"^). In other words, an efficient "boost" of Lya 
with respect to the continuum would require such extreme 
ISM conditions. 

A high dust content has to be embedded in 
clumps: As explained in Sect 3.2.2 and shown in Fig. 
18, no enhancement of EW(LyQ;) is found below a certain 
critical dust content (noted Tc in Fig. 6). A high dust con- 
tent is indeed needed in the ISM in order to absorb more 
efRciently UV continuum photons than Lya photons, which 
thus produces an enhancement of EW(Lya). As shown in 
Fig. 18, a colour excess higher than E{B — V)caizetu = 0.32 
(ie. Tc ~ 5) would be needed in order to enhance EW(Lya) 
by a factor higher than 3. This dust content limit mainly 
depends on the covering factor CF of the ISM, where it 
decreases as CF decreases. 

The distribution of the clumps around the 
stars: The spatial distribution of the clumps around 
the stars plays an important role in the enhancement of 
EW(Lya). A covering factor CF close to unity is needed in 
order to get a enhancement of EW(Lya) as high as those 
shown in Fig. 18. The enhancement of EW(LyQ;) is indeed 
maximized when CF — 1, but it strongly decreases as 
CF decreases around the stars. A low covering Factor CF 
does not allow to block effectively UV continuum photons 
than Lya photons, thereby decreasing the enhancement 
of EW(Lya). In particular, under the physical conditions 
adopted in Fig. 18, we notice that the enhancement of 
EW(Lya) stays lower than 1.38 for CF <, 0.68. 

Let's mention that, in both figures 18 and 19, the to- 
tal Lya flux is taken into account to derive the EW(Lya) 



enhancement. However, an observer could measure a lower 
EW(Lya) in practice if Lya is scattered into an extended 
low surface brightness region, as recently observed around 
distant starburst galaxies (Steidel et al. 2011) or in nearby 
galaxies (Ostlin et al. 2009). 

Concerning the required dust extinction and the high 
density contrast of neutral hydrogen (i.e. njc/nc < 10"''), 
these conditions seem more characteristic of molecular 
clouds embedded in a very hot and ionized medium. The 
Neufeld model could therefore only work in such a galac- 
tic environment. It is nevertheless interesting to notice 
that these necessary conditions to get an enhancement of 
EW(Lya) in a clumpy ISM are in perfect agreement with 
those originally predicted by Neufeld (1991). In his original 
paper, Neufeld proposed indeed three suitable conditions 
for an enhancement of EW(Lya) in a clumpy ISM: the 
interclump medium (ICM) must exhibit a very low den- 
sity with a negligible small absorption and scattering coef- 
ficients to Lya photons, the clumpy ISM must show a large 
covering factor (CF) and a sufficiently small volume filling 
factor that the ICM can "percolate" (i.e. every part of the 
interclump medium must be connected to every other part) , 
and the probability for Lya photons of being reflected by 
the clumps must be higher than the probability to be ab- 
sorbed. In particular, this last probability should increase 
when Vcxp decreases and iVni increases, which is well consis- 
tent with our results. In addition to the criteria proposed by 
Neufeld (1991), our work highlights the extreme sensitivity 
of the Neufeld scenario to the kinematic of the clumps (i.e. 
the large scale outflows and the velocity dispersion of the 
clumps), and its strong dependence on the dust content. 
Given the ubiquitous evidence for outflows from most star- 
forming galaxies and the widespread presence of dust, it 
is essential to take these effects into account. Furthermore, 
beyond the work of Neufeld (1991), our detailed radiation 
transfer models including the Lya line but also transfer of 
continuum photons at other wavelengths, also allow us to 
predict consistently the resulting attenuation (reddening) 
and the detailed Lya line profile, which can directly be 
compared to observations. 

5.3. Studying the ISM through the Lya line profile 

Although the Lya line profiles emerging both from homoge- 
neous and clumpy ISMs are quite similar (Fig. 13), we have 
identified two main effects produced by the ISM dumpi- 
ness on the Lya line profiles. First, we have seen that an 
extremely clumpy ISM favours the formation of a peak at 
the line center of the line profile (wobs=0 km s~^). Second, 
since Lya escape is facilitated in a clumpy medium, the ef- 
fect of the dust on the Lya line profile is less efficient than 
in homogeneous media (Figs. 15 and 16). This can lead to 
intense Lya emission emerging from a very dusty clumpy 
ISM, whereas an absorption line profiles would emerge from 
a homogeneous medium with the same dust content. 

This second effect of the ISM dumpiness on the Lya line 
profiles can be a source of uncertainties if we aim to derive 
some informations on the ISM of distant starburst galax- 
ies (kinematics Vcxp, Hi column density TVhIj dust content 
Ta) from Lya line fitting (Verhamme et al. 2008; Schaerer 
& Verhamme 2008). Indeed, depending on the homogene- 
ity/dumpiness of the ISM, different derived parameters of 
the ISM can be obtained studying the same sample of galax- 
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ies. Among the ISM parameters it is possible to derive from 
the fit of both the Lya line and the UV continuum (voxp, 
■^Hi, T^), the dust optical depth (ja) seems the most un- 
certain parameter, given the dependence on the degree of 
clumping. On the positive side, the expansion velocity may 
probably be best determined. In asymmetric Lya profiles 
the frequency of the second bump - if present - traces quite 
well 2 X Vcxp, both in homogeneous or clumpy media, as al- 
ready pointed out by Verhamme et al. (2006). Hovifever, 
complications may be that two bumps (i.e. those located at 
the line center and those redshifted at 2 x Vcxp in Fig. 16) 
are not detectable (e.g. due to insufficient spectral resolu- 
tion), and the second bump (i.e. shifted at 2 x Vexp) may 
not always correspond to the peak of the profile. Also, de- 
pending on the column density the distinction of the two 
peaks may not be very easy (Verhamme et al. 2006). 

In any case, compared to the typical Lya line profiles 
observed in distant Lya emitters and LBGs, we note that 
most of them are asymmetric lines with redshifted peaks, 
which seem difficult to reconcile with the profiles predicted 
for very clumpy shell geometries, both static or expand- 
ing, since the peak is then expected to show negligible red- 
shift. Furthermore, Lya absorption lines are usually found 
among the reddest (presumably more dusty) LBGs, and 
the Lya equivalent width correlates with reddening (e.g. 
Shapley et al. 2003) , facts naturally explained by radiation 
transfer models with a homogenous ISM (Verhamme et al. 
2008; Schaerer & Verhamme 2008). These findings also ar- 
gue against a very clumpy, high-contrast ISM, at least for 
the majority of LBGs. 

5.4. Comparison with Hansen & Oh (2006) 

We now compare our results with the recent numerical 
study of Lya transfer in multiphase and dusty media from 
Hansen & Oh (2006) (hereafter HO06). The clumpy me- 
dia studied in HO06 are only extremely clumpy, that is 
composed of very optically thick spherical clumps (Hi + 
dust) distributed within an empty interclump medium. 
From such clumpy media, HO06 deduce analytical formu- 
lae fitting the behavior of both the continuum and the Lya 
escape fractions as a function of gas geometry, motion and 
dust content. In summary, assuming an isotropic scattering 
by dust grains {g — 0) in optically thick spherical clumps, 
the fitting formula is a function of two parameters (Eq. 59 
in HO06): 

fescM = l/cosh(V2ec(j^)iVo) (15) 

where Nq (a geometrical parameter) corresponds to the av- 
erage number of clumps encountered by photons before es- 
caping the medium in the absence of absorption, and where 
ec(i^) (a dust parameter) corresponds to the probability of a 
photon of frequency v to be absorbed rather than reflected 
by a clump. In the clumpy shell geometries of our model, 
both parameters can be derived in the following way. Firstly 
we notice that Nq tends to evolve as Nq = fc^ + 1.06/c, 
where fc corresponds to the mean number of clumps inter- 
sected along a random line of sight. The parameter edi') 
can simply be estimated using the formula derived in HO06 
(eq. 27 of their paper): 

ee(;.) = 2^)/{l + v/^) (16) 

with e(i^) as the absorption probability per interaction (HI 
or dust) at frequency v. In particular, e{i>) is thus given by 
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Fig. 20. Test of the equation 15 of HO06. We show here the evo- 
lution of the UV escape fraction from two different extremely 
clumpy shell geometries as a function of the absorption param- 
eter ec. The first medium (square dots) shows No = 1.00 and is 
buih with FF = 0.10, CF = 0.38 and mc/nc = 0. The second 
medium (circle dots) shows A'^o ~ 4.28 and is built with FF = 
0.15, CF — 0.80 and nic/nc ~ 0. In both media, each clump 
is opaque to dust extinction. The clumps are therefore optically 
thick for UV continuum photons. Both full lines represent re- 
spectively the equation 15 for A'^o = 1.00 (top line) and A^'o = 
4.28 (bottom line). 

e(i/) = 1 — a for the UV continuum photons, with a the 
dust albedo. 

We compare in Fig. 20 the equation 15 (full lines) 
to the UV continuum escape fraction derived from the 
clumpy shell geometries of our model (grey dots). In this 
figure, two clumpy shell geometries are studied. The first 
structure, buih with FF = 0.10, CF = 0.38 and nic/nc 
— 0, shows A^'o ~ 1.0 (top line). The second structure, 
assuming FF — 0.15, CF = 0.80 and nic/nc — 0, shows 
A'o = 4.28 (bottom line). In Fig. 20, we just change the 
values of edi^) changing the albedo a of the dust grains, as 
shown in Eq. 16. Although a certain difference appears as 
Cc tends towards 1, the continuum escape fraction deduced 
from our numerical simulations are rather well fitted by 
the analytical formula of HO06 (Eq. 15). In particular, the 
UV escape fraction follows well the same dependance on 
A^o and Sc as predicted by Eq. 15 for the low absorption 
regime (e^, < 0.4). The difference observed close to = 1 
is also observed in HO06 and is explained by a geometrical 
effect. As suggested by these authors, a better fit can be 
obtained in this regime rescaling the term ec{v)NQ (Eq. 
15) as k€c{v)Nq^ where k is a unity fitting parameter. 

Although our numerical simulations reproduce reason- 
ably well all the results and fitting formula of HO06, 
our conclusions diverge from theirs concerning the Neufeld 
model. As an application of their numerical study, HO06 
give a quantitative estimation of the EW(Lya) enhance- 
ments produced by different clumpy and dusty ISMs (Fig. 
18 of their paper). We reproduce in Fig. 21 (left panel) the 
results of HO06. These results are based on the following 
assumptions concerning the clumpy media: 1) the clumps 
are extremely opaque to Lya photons, but 2) each clump 
is not opaque to dust extinction (i.e. the dust extinction 
optical depth for a single clump is Tq < 1) and are there- 
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Fig. 21. Evolution of the Lya line-to-continuum escape fractions ratio, /esc(LyQ!)//csc(UV), as a function of the dust optical depth 
Ta in six different extremely clumpy shell geometries. The left panel illustrates the results of HO06 (Fig. 18 of their paper), whreas 
the right panel shows those of our own simulations. Note the different vertical scales. Each clumpy shell geometries are characterised 
by different values of A^o and are built with different set of parameters. Both structures showing A^'o = 1 (dotted lines) are built 
adopting FF=0.10, CF=0.38, nic/nc=0, T=10** K, iVm=10^^/cm"^ and Vexp=50, 250 km s"^ Both structures showing A^o = 4 
(dashed-dotted lines) are built adopting FF=0.15, CF=0.80, nic/nc=0, T=10'* K, ]Vm=10^^ cm"^ and v^xp = 50, 250 km s~\ 
Finally, both geometries showing A''o = 10 (sohd lines) are built with FF=0.13, CF=0.90, nic/nc=0, T=10* K, iVm = 10^^ cm~^ 
and Vexp = 50, 250 km s~^. In the left panel, each curves are obtained assuming /csc(UV) = e"'^" and, respectively, /esc(Lya) — 
(0.94; 0.68; 0.38) for 50 km s"^ and (0.81; 0.43; 0.18) for 250 km s"^ Conversely, all curves shown in the right panel are derived 
from full Monte Carlo simulations. Like in HO06, the Lya escape fraction is derived in our simulations at the line center (ie. the 
intrinsic Lya line is described as a delta function in all simulations). 



fore optically thin for UV continuum photons. Given these 
assumptions, HO06 adopt a constant value of /osc(Lya) 
(which is derived by using their analytical approximations 
based on equation 15), whereas the UV continuum escape 
fraction is assumed to behave Hke in homogeneous media. 
In a homogenous medium composed of dust grains with 
an absorption albedo « 0.5, the UV escape fraction is 
approximately given by the following equation: 

/csc(UV) = 1/ cosh(V4ed(T2+T,)) (17) 

On the right panel of the Fig. 21 we show the results 
obtained studying the same media than HOGG, but us- 
ing our own numerical approach. In particular, we derive 
each curves of the right panel studying the real evolution 
of /esc(UV) and /esc (Lya) with full Monte Carlo simula- 
tions. We can clearly see that we cannot reproduce, neither 
quantitatively nor qualitatively, the curves of HOGG. This 
strong difference is explained by both assumptions made by 
HO06 on /csc(UV) (i.e. Eq. 17) and /esc (Lya) (i.e. a con- 
stant value under any dust optical depth t^) , which do not 
reproduce the correct evolution of /esc (Lya) and /osc(UV) 
in extremely clumpy media. We compare in Fig. 22 the UV 
and Lya escape fractions deduced from our simulations to 
the assumptions made in HO06. These assumptions clearly 
underestimate the correct values of /csc(UV) and overesti- 
mate those of /esc (Lya). This explains why the enhance- 
ments of EW(Lya) obtained in HOGG (left panel in Fig. 
21) are much higher than those deduced from our simu- 
lations (right panel). Furthermore, the right panel shows 
the inversion of the curves (in terms of iVp) as increases. 
The highest enhancement of EW(Lya) are indeed produced 
by clumpy media containing the highest number of clumps 
around stars. 



Drawing the parallel with our study of the Neufeld sce- 
nario in Sect. 5.2, we can notice that most of the mod- 
els showing an enhancement of EW(Lya) in Fig. 21 (right 
panel) respect quite well each conditions under which the 
Neufeld scenario works in a clumpy ISM : 1) the interclump 
medium is optically thin for Lya photons {nic/nc — 0), 2) 
the expansion velocities Vexp is slow and uniform and 3) 
the enhancements of EW(Lya) occur above a certain dust 
optical depth Tc equal to (0, 0, 0.3) in each clumpy me- 
dia (A'o = 1,4, 10). Furthermore, as expected when an en- 
hancement of EW(Lya) occurs in a clumpy ISM, the Lya 
line profiles emerging each clumpy geometries studied in 
Fig. 21 are symmetric and peaked at the line center, as 
expected. 

5.5. Recent models of Lya transfer in clumpy large-scale 
outflows 

Recently Dijkstra & Kramer (2012) have presented ra- 
diative transfer calculations of Lya photons propagating 
through clumpy, dusty, large scale outflows using phe- 
nomenologically motivated models constrained by absorp- 
tion line measurement from Steidel et al. (2011). The cal- 
culations of Dijkstra & Kramer mostly focus on the Lya 
surface brightness distribution. However, since their calcu- 
lations do not follow the behavior of the UV continuum, it is 
not possible to infer any information on the presence or ab- 
sence of an efficient "Neufeld effect" . Furthermore the Lya 
line profiles predicted by their models are not presented. 
It is therefore not possible to compare our results with the 
calculations of Dijkstra & Kramer (2012), and to confront 
their model to the most direct and sensitive observable, the 
Lya line profile itself. 
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Fig. 22. Predicted UV and Lyo? escape fraction as a function 
of Ta for the clumpy media studied in Fig. 21 (A^o = 1 (trian- 
gle), iVo = 4 (circle) and iVo = 10 (square)). Top: /csc(UV) ob- 
tained from our simulations compared to the equation 17 used 
by HO06 to carry out their figure shown on the left panel of 
Fig. 21. Middle: Lya escape fraction /csc(LyQ!) predicted by our 
models for Vexp = 50 km s^^. Bottom: Same as middle panel 
for Vexp = 250 km . Like in Hansen & Oh (2006), the Lya 
escape fraction is derived here at the line center (ie. the intrinsic 
Lya line is a delta function in our simulations) While a clear de- 
crease of the /esc (Lya) with increasing is found, HO06 rather 
assumes a constant value of /esc (Lya), indicated in the figure, in- 
dependently of the dust optical depth t^. More precisely, HO06 
assumes /esc(Lya) = (0.94; 0.68; 0.38) for Vexp ~ 50 km s~^ and 
/esc(Lya) = (0.81; 0.43; 0.18) for Vexp = 250 km s-^ 



6. Summary and conclusions 

To examine and understand the effects of clumpy ISM 
structures on the Lya line and UV observations of star- 
forming galaxies we have carried out detailed radiation 
transfer calculations using our 3D Lya Monte Carlo code 
MCLya (Verhamme et al. 2006; Schaerer et aL 2011). 
Indeed, clumping can in principle significantly alter the 
transfer of Lya in galaxies, as shown early by Neufeld 
(1991), and has often been invoked to explain strong Lya 
emission or a higher transmission for Lya photons than for 
the UV continuum (e.g. Kudritzki et al. 2000, Malhotra & 
Rhoads 2002, Rhoads et al. 2003, Shimazaku et al. 2006). 
However, only few detailed numerical studies of these ef- 
fects has so far been undertaken (Haiman & Spaans 1999; 
Richling 2003; Hansen k Oh 2006; Laursen et al. 2012)^ 
albeit with some simplifying assumptions in most of these 
works. Furthermore we wish to identify in which conditions 
clumping affects the line transfer and how much, and how 
this is reflected in the emergent Lya line profiles. 

Our radiation transfer calculations allow us to study 
simultaneously the dependence of both the Lya and the 
continuum escape fractions, the Lya equivalent width 



EW(Lya), and the Lya line profiles on the Hi content, 
the dust content, kinematics, gas geometry, and clumping 
properties. Since spherically symmetric outflows with a ho- 
mogeneous Hi shell are able to reproduce a large variety 
of observed Lya line profiles in Lyman break galaxies and 
Lya emitters (Verhamme et al. 2008; Schaerer & Verhamme 
2008; Dessauges-Zavadsky et al. 2010), the same geometry 
is used to study how a clumpy ISM structure alters the 
Lya line and UV continuum. This clumpy geometry is also 
chosen since it has been shown to reproduce observable con- 
tinuum properties of starburst galaxies and the Calzetti at- 
tenuation law (Gordon et al. 1997; Witt & Gordon 2000; 
Vijh et al. 2003). 

Our main results can be summarized as follows: 

— A clumpy and dusty medium is always more transparent 
to Lya, UV and optical continuum photons compared 
to an equivalent homogeneous medium of equal dust 
content, as already known from earlier studies (Boisse 
1990; Hobson & Scheuer 1993; Witt & Gordon 1996, 
2000). A clumpy medium thus allows to decrease the 
global effect of the dust absorption on any radiation. 

— The UV and optical continuum escape fraction depend 
on three parameters in homogeneous and clumpy shell 
geometries: the dust content, the "dumpiness" of the 
dust distribution (described by the density contrast 
nic/nc), and the covering factor CP (defined here as 
the fraction of solid angle of the central photons source 
covered by the clumps). In a general way, the continuum 
escape fraction decreases as the dumpiness of the dust 
distribution decreases (i.e. njc/nc increases), and as 
both the dust content and the covering factor increase. 

— Three additional parameters, i.e. in total six, control the 
Lya line transfer: the dust content, the density contrast 
nic/nc, the covering factor, as well as the Hi column 
density, the velocity field and the gas temperature. The 
Lya escape fraction always increases with increasing the 
dumpiness (i.e. njc/nc decreases), and with decreas- 
ing dust content and covering factor. However, the Hi 
column density and the kinematics of the gas do not 
affect the Lya escape in the same way for homogeneous 
or clumpy media. That creates two different regimes for 
the Lya radiative transfer in clumpy media. 

— The first regime (called "low contrast" regime in this pa- 
per) comprises homogeneous and weakly clumpy shell 
geometries, corresponding to an interclump density 
above nic/riQ ^ 1.5 x 10^'' for the physical condi- 
tions adopted in our model (such that the interclump 
medium is also optically thick for Lya). In this regime, 
the Lya escape fraction increases with increasing ex- 
pansion velocity Vexp and with decreasing Hi column 
density A'hi, as for a homogenous ISM. The Lya escape 
fraction is then always less or equal to the UV escape 
fraction, which implies that the emergent Lya equiva- 
lent width i?VFobs(Lyck) is lower than the intrinsic one 
EW,^t(Lya). 

— The second regime (called "high contrast" regime) is 
found in the most extremely clumpy shell geometries of 
our model {nic/nc ^ 1.5 x 10~^ for the physical con- 
ditions adopted in our model). This corresponds to a 
clumpy medium composed of very dense clumps embed- 
ded in an interclump region which is optically thin for 
Lya photons. Two main differences appear compared to 
the other regime. First, as was originally suggested by 
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Neufeld (1991). it is possible to observe a Lya escape 
fraction which is higher than for the UV continuum. 
In particular, this is possible above a certain "critical" 
dust optical depth Tc- Second, whereas for Tq < Tc the 
Lya escape fraction behaves as in the "low contrast" 
regime, the opposite behavior is found for high enough 
dust content {j^ > Tc), where /esc (Lya) increases with 
decreasing Vcxp and inereasing Nm. 

— Overall wc have identified two main effects of the ISM 
dumpiness on the shape of the Lya line profiles. First, 
extremely an clumpy ISM favours the formation of a 
peak at the center of the line profile (wobs=0)- Second, 
the intensity of the Lya line increases as the dumpiness 
of the medium increases, as expected. 

— Schematically, the following Lya line profile mor- 
phologies are predicted from homogeneous and clumpy 
shell geometries: "double-peak" profiles with identi- 
cal/similar peaks symmetric around the source redshift 
(for static media), asymmetric redshifted profiles (from 
expanding media), and absorption line profiles (from 
very dusty, homogeneous or weakly dumpy media). 
These types have already been identified in homoge- 
nous models (Verhammc et al. 2006). In very clumpy, 
static and dusty shells, a new category is found: "three 
peaks" profiles similar to the double-peak profiles with 
an additional third component at line center. 

As an application of our study, we have examined the 
conditions under which the Neufeld model (Neufeld 1991) 

can work in a chimpy ISM, i.e. when an enhancement of the 
observed EWobs(Lya) can be obtained. We find that the 
following five conditions must be simultaneously fulfilled 
for the "Neufeld" effect to be effective. 

— The Hi content must be very low between 
clumps, typically the Hi column density in the in- 
terclump region must be 3 x 10^"* cm~^ to remain 
optically thin for Lya photons. Otherwise, Lya pho- 
tons scatter strongly against Hi atoms localized between 
clumps and cannot escape the clumpy medium more 
easily than UV photons. 

— The galactic outflow has to be slow with outflow 
velocities of the order of Vexp 2 x FW H Mint{Lyct) 
km s~^. Otherwise Lya photons are too redshifted from 
the chimps and cannot scatter anymore on the surface 
of clumps, as suggested in the Neufeld model. 

— The galactic outflow has to be as uniform and 
constant as possible in velocity for efficient inter- 
actions of Lya photons with dust. 

— A high dust content must be embedded in 
clumps to absorbe as much as possible the UV con- 
tinuum, which increases the Lya equivalent width. For 
the physical conditions and clumpy shell geometries 
adopted here, we find that an enhancement of EW(Lya) 
by a factor 3-4 can only occur for a colour excess E(B- 
V) Z 0.3. 

— A large covering factor is needed in order to get a 
noticeable enhancement of EW(Lya). 

The above conditions are in agreement with the general 
findings of Neufeld (1991). However, our results differ from 
those of Hansen & Oh (2006), who make some simplifying 
assumptions, which are not consistent with more rigorous 
radiation transfer calculations. In our study, the Neufeld 



model does not work as easily as suggested in Hansen & 
Oh (2006). 

Given our results it seems quite unlikely /difficult to 

find conditions in a clumpy, spherically symmetric ISM, 
where Lya photons can escape more easily than the nearby 
UV continuum, i.e. where the phenomenon suggested by 
Neufeld (1991) can be at play. Furthermore, when these 
conditions are fulfilled we generally find that the emergent 
Lya line profile shows emission at line center and little 
asymmetry. Such profiles do, however, not represent the 
profiles typically observed in high-redshift galaxies, which 
are known to be redshifted in the galaxy rest-fame and 
asymmetric. Other arguments against the Neufeld effect 
being effective may be if the sites of Lya emission arc rel- 
atively close to or within cold, dusty environments such 
as molecular clouds, which could absorb more efficiently 
Lya photons than the UV continuum (Laursen et al. 2012; 
Verhammc et al. 2012). 

The simulations from this paper and the success of ho- 
mogenous, spherically expanding models in reproducing the 
large variety of observed Lya line profiles and velocity shifts 
between photospheric, low ionization absorption, and the 
Lya line (Verhammc et al. 2008; Schaerer & Verhamme 
2008) seem to indicate that effects to due inhomogeneities 
in the ISM and deviations from spherical geometry are not 
dominant. Why simple geometries work so well may ap- 
pear somewhat puzzling, and certainly remains worth un- 
derstanding more in depth. More detailed observations and 
sophisticated radiation transfer models within multi-phase 
ISM, high-resolution simulations may help to shed more 
light on this question and to examine how well the physical 
parameters derived from Lya line profile fits using simple 
spherically expanding model represent reality. 
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